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Abstract 

We propose a free field representation for tlie form factors of descendant operators in tlie BuUougli- 
Dodd model. This construction is a particular modification of Lukyanov's technique for solving the 
form factors axioms. We prove that the number of proposed solutions in each level subspace of the 
chiral sectors coincide with the number of the corresponding descendant operators in the Lagrangian 
formalism. We check that these form factors possess the cluster factorization property. Besides, we 
propose an alternative free field representation which allows us to study analytic properties of the 
form factors effectively. In particular, we prove that the form factors satisfy non trivial identities 
known as the "reflection relations". We show the existence of the reflection invariant basis in the 
level subspaces for a generic values of the parameters. 

o ; 1. Introduction 

<~i . The form factors of local operators in two dimensional integrable models can be exactly obtained as the 

solutions to a system of functional equations known as the form factor axioms PQ [5] . These equations 
specify the correct analytical properties of the form factors predefined by the scattering matrix of the 
model. In principle, one can solve these equations recursively. The space of solutions to the form factors 
axioms is supposed to give a faithful representation of the operator content of the model. However, the 
correct identification of the form factors of the specific operator within tliis space of solutions is in general 
non trivial problem. In the case of the exponential operators there exist two useful criteria which assist 
in the identification. The first one restricts the asymptotic behavior of the form factors at large values 
of rapidity [4J while the second one is the so called cluster factorization property fSj. But the family of 
exponential operators are far from exhausting the full set of operators in the theory, which also contains 
descendant operators. The two criteria we discussed are not sufficient to solve the identification problem 
for these operators. Therefore, we need an effective approach to study the form factors of descendant 
operators. 

Here we consider the Bullough-Dodd model [6l [7] . A convenient method for solving the form factor 
axioms is the free field representation proposed by Lukyanov [5]. By using this method the solutions to 
the form factor axioms in the Bullough-Dodd model that correspond to the exponential operators were 
found and completely identified [H]- Following the guidelines of [TD] we find the free field representation 
for the form factors of descendant operators. However, the problem of identification is not solved in this 
case. 

The proposed free field representation allows us to study analytic properties of the form factors 
efficiently. In particular, we prove that the form factors satisfy the reflection relations [llj. These 
relations originate from the similar relations in the Liouville theory and establish connection between 
the exponential operators e""^ with different values of the parameter a [12]. Besides, we prove the 
existence of the reflection invariant bases in the Fock spaces. 

This paper is organized as follows. In Section [2] we describe the Bullough-Dodd model and specify 
its operator contents. In Section [3] we briefly describe Lukyanov's free field representation for the form 
factors of exponential operators. We introduce auxiliary construction and find the free field representation 
for the set of solutions to the form factor axioms. The number of these solutions on each level subspace 
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of the Fock module is shown to be coincide with the dimension of the corresponding level subspace in 
the Lagrangian formalism. We prove that these solutions satisfy the cluster factorization property and 
study their asymptotic behavior. In Section 0] we introduce alternative free field representation which 
is used to obtain recurrence relations between the form factors and prove the reflection properties. In 
Section [5] we consider an explicit example, namely, level (2,0) descendant operators. Finally, we draw 
our conclusions in Section [51 



2. The Bullough Dodd model 

The Bullough-Dodd model is defined by the Eucledean action [6 , 7] 



Sbd — 



d^x (^j^idu^f + /^(e^*"" + 2e-Vf '^)^ . (2.1) 



Here b is the coupling constant and n is the regularized mass parameter. The spectrum of the model 
consist of the single particle A of mass M which can be expressed in terms of the parameter fi ^] . This 
particle appears as the bound state of itself in the scattering processes 

AxA^A^AxA. (2.2) 

Hereinafter we shall use the following notation 

Q^b + b-\ Lu^e'''/^. (2.3) 

Also, we shall use the light-cone coordinates 

1 - liOa^ 5^ 
z~x — X , z — x +x , o——, o — —. 

oz az 

The integrability of the model implies that the scattering processes are purely elastic [13]. Therefore, 
the n-particle ^-matrix factorizes into the 71(71 — l)/2 two-particle scattering amplitudes |14| 

tanh i(6i + 2hl) tanh i(6' - 1^) tanh i(6i - 2i^) 

SiO) = ^- ^ ^- 1^^^ ^- (2.4) 

^' tanhi(6i- 2HL)tanhi(e'+ |^)tanhi(6i+ ^) ^ ' 

The ^-matrix is invariant under the weak-strong coupling constant duality b — > b^^ . A simple pole of 
the scattering matrix located at 6' = 27ri /3 corresponds to the bound state of the particle A itself in the 
scattering process (j2.2l) . 

Let us consider the operator contents of the model. The space of local operators of the Bullogh-Dodd 
model consists of the exponential operators 

Va{x) = e"'^^^) (2.5) 
and their descendants, i.e. the linear combinations of the fields 

5"V • ■ • . . . 9"^(^e"^(^\ (2.6) 

Here the pair of integers (n, n) given by 

r s 

n = ^^ni, n — '^^nj (2-7) 

*=i 3=1 

is called the level of descendant operator. The numbers n and n separately are called chiral and antichiral 
level correspondingly. Any exponential operator is characterized by its scaling dimension Aq in the 
ultraviolet regime which is related to the parameter a. The scaling dimension of the corresponding 
descendant operator at the level (n, n) is given by -I- n -I- n, while its spin is 71 — n. 

Let us consider the radial quantization picture at some point in the Eucledean plane, e.g. a; = 0. In 
the vicinity of this point the field Lp{x) can be expanded as 

^{x) = Q - log zz + y ^z-'" + y ^z-™. (2.8) 
^ im ^ im 
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Here the operators Q, P, a,„ and a„j form a Heisenberg algebra with the fohowing commutation relations 
[P, Q]=-i, [a„i,a„] = m5„i+„^o, [a„i, a„] = m(5„i+„^o- (2.9) 

In the radial quantization picture the exponential operator Vq(0) corresponds to the highest weight vector 
\a)rad defined by the relations 

ahn\a)rad = am\a)rad = (m > 0), P\a)rad = a\a)rad, \a)rad = e''^\vac)rad (2-10) 

Let J-a be the Fock module spanned on the vectors generated by the elements a_m (m < 0) acting on 
the highest weight vector \a)rad- Similarly, the is a Fock module spanned on the vectors generated 
by the elements a_„i (m < 0) acting on the same highest weight vector. Evidently, the modules and 
T are isomorphic. The space of states is given by the tensor product J^a (X) of chiral and antichiral 
components. Up to some numerical factor the descendant operators (|2.6p correspond to the vectors 

a_„^ . . . a_„^a_sj . . . a_s^ \a)rad (0 < ni < . . . < n^, < fii < . . . < fis). (2.11) 

Each module Ta admits a natural gradation into the level subspaces J-a = ®J^o-^a,"j where each level 
subspace J-a.n is spanned on the vectors (|2.1ip with h — and J2 = The dimensions of these 
subspaces are given by the following generating function 



^g"dim^,„=nT^- (2.12) 

n=0 m=l ^ 



The dimension of the level (n, n) subspace if given by the product dim J^a,„ • dim Ta,n- 

3. Free Field representation for form factors 
3.1. Form factors of exponential operators 

Let us briefly describe Lukyanov's free field representation for the form factors of exponential operators 
in the Bullough-Dodd model [9\. First, we consider a pair of operators A+(6'), A^(6') and define the one 
and two-point trace functions as follows 

((A-(0)))=1, {{A'^'{0')A'^{0))) = [R{0-eT'\ 'T',a = ±, (3.1) 
where R{6) is the two-point minimal form-factor given by |15| 

/ P°° df cosh I sinh sinh jQf \ 

i?(0) exp(-4 / f ^__!Q^cosh(t-^) . (3.2) 

^ Jo t smh t smh ^ tt / 

The multi-point functions can be calculated by means of the Wick's theorem. Let us define the normal 
ordering procedure : . . . : by the following relation 

A^"(0w)...A"n^i)=:A-«(0^)...A-H^i): H ((A"^' (^j)A"' (^O))- (3.3) 

l<i<j<N 

Hereinafter we shall use the notation 

K\e)^:K+{e--)K-{e +'-):. (3.4) 

The Lukyanov's generators are given by 

T(6') =p(e"PA+(6i) + e-''^PA-(6i) + /iA°(6i)), (3.5) 



where 



2s n-^ ^, 

eg ' 

6Q 2' 

"°° dt cosh I sinh ^ sinh 3^ X 



(3.6) 



/ sinf f f 

\ sm^sm^ V Jo 



l5iii ^ J ^ sinli t cosli 2 
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The form factors of exponential operators are given by the muhi-point trace functions, namely 

(mc|K(O)|0i,...,f?jv) = {e''nfa{Oi,...,eN) = (e'^'^)((r(^jv) • • (3-7) 

where (e"'^) is the vacuum expectation value of the corresponding exponential operator [llj. The 
functions /q(0i, . . . ,0n) are analytic functions in the variables 9i with complicated analytic structure. 
Taking into account p.ip we immediately conclude that these functions can be represented as follows 

Uei,...,dN)^p''JN,a{e'' ■■■,e'") H R{e,-e/). (3.8) 

i<i<j<N 

Here the functions JN,a{xi, . . . ,xn) are symmetric rational functions in the variables Xi = Their pole 
structure is governed by the form factor axioms. The only poles located at relative rapidity difference 
6ij — in and 9ij = 27rz/3 will be referred to as the kinematical and the bound state poles correspondingly. 



3.2. Form factors of descendant operators 

We argue that the form factors of descendant operators possess the free field representation that is very 
similar to Lukyanov's one p.7p . The main idea is to modify the generators T{6). Let us consider a 
commutative algebra A = Q^^qAu generated by the elements {a_„} with n > 0. Each level subspace 
An is spanned on the elements Y[ ct-ui such that '^Ui = n. We also introduce another copy A of the 
algebra A generated by the elements {q;_„}. The canonical homomorphism between these algebras is 
defined by the following relation: for any h ^ A\et us define h ^ A according to the rule q;_„ -t> Oi-n- 
The element g = hh' will be referred to as the level [n, n) descendant ii h £ An and h' E An- 
Let us define a bracket on the algebra A by 

/ oc oo \ oo 

n n = n kjA^,^ (3.9) 



and consider the currents 



\7n — l m—1 



a{z) — exp [ a-r 

m—l 

00 



5(z)=exp( - }_^a_„(-zrl, (3.10) 

/ oo 

c(z) = exp f ^ {oj-'" - {-iruj"')a.,nz 

^ m=l 

By using these currents we modify the Lukyanov's generators as follows 

= p(e*"f a(e^)6(e-'')A+(0) + e-'^'P b{e^)aie^'')A~ (0) + hc{e^)cie~^)A"{e)) . (3.11) 

For any element g E A® A consider the function 

/f (01, ...,9n)^ (meN), . . . , m))), 5). (3.12) 

It is straightforward to check that this function is a solution to the form factor axioms. Indeed, the 
Watson's theorem and the crossing symmetry condition are evidently satisfied while the kinematical 
and the bound state pole conditions are satisfied if the currents (I3.10p are subjected to the following 
conditions 

a{z)b{—z) = 1, c(z) = a{zuj^^)b{zuj), 

which are the case. Consequently, the function (j3.12l) determines the form factors of an operator from the 
Fock space (J^ ® T)Va{x). The descendant operator which corresponds to this function will be denoted 
by v/{x) that is 

{vac\V,^ix)\ei, . . . , 0^) = {e'^nf^iOi, ■ • ■,0n). (3.13) 
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From expression p.l2p it follows that the solutions to the form factor axioms are given by 

N 

/f (01, . . . , ^at) = J^, Je^S . . . , e^«) n m ~ Oj). (3.14) 

i<j 

Here the functions J^{xi, . . . , xn) are symmetric rational functions in the variables x,; = e^' with 
the kinematical and bound state poles located at rapidity differences 9ij = in and 9ij — ±2i7r/3 
correspondingly. Evidently, the J functions which correspond to the exponential operators are given 
by "^N a ~ JN,a- It IS Straightforward to get an explicit expression for the function ^, namely 

J%^^{x,,...xm)= /i^'°e(#^+-#^-)'^fPnX+|X_|Xo)x 

/++/_+/o=/ 



where 




In expression (j3.15|) we introduce a set of integers / = {1,...,A^} and the sum is taken over all 
decompositions of the set / into the three subsets I„, a = {+,—,0}, such that /+ U U /q — I 
and /(j/ n /ct = if cr' 7^ 0". With every subset I^r we associate the subset Xa- = {xi\i G la}- The 
functions P^{X\Y\Z) are polynomials defined by the following relations 

P"--{X\Y\Z) = S,n{X) - {-irS,n{Y) + (c^-" - (-l)"c^")5„(Z), 

P^-"^{X\Y\Z) = S.,n{Y) - {-irS-,n{X) + {cj-"' - (-1)"W")5_„(Z), (3.17) 

p9i92 ^ pgi pg2 ^ pcigi+c2a2 ^ ^^pai ^^pg2 (V51, c,2 G «) A ci, G C). 

Here we denote the power sums of order m by 

N 

5,„(xi,...,xjv) = (3.18) 

1=1 

In summary, in the framework of the free field representation we obtain the set of solutions to the 
form factor axioms. Besides, in (I3.14p and (I3.15P we give explicit expressions for these functions. The 
form factors of exponential operators which correspond to the case g — 1 have been studied in detail 
We want to prove that the space of proposed solutions can be bijectively mapped into the space of 
descendant operators in the Lagrangian formalism. In other words, we want to prove that for any 
descendant operator in the Lagrangian formalism there exist a solution to the form factor axioms given 
by certain linear combination of the functions 

It is a challenging problem to make an identification between the states from the space A® A and 
those ones from T ® T. Certain requirements are necessary to identify the form factors of a specific 
operator among all these functions. In the case of exponential operators the solution to this problem is 
known. It is sufficient to impose restriction on the asymptotic behavior of the form factors at large values 
of rapidities [3] and demand the cluster factorization property to be satisfied [5j. However, additional 
criteria are required in the case of descendant operators. This problem is not solved in general cases. 
Nevertheless, let us consider the analytic properties of the proposed set of functions. 

3. 3. Cluster factorization property and asymptotic behavior 

We consider the element g = hh' such that h G A, h' G A. We are interested in the asymptotic behavior 
of the function which corresponds to this element, namely 

ffiOi, 0„, 6.^+1 +A,9n+ A), (3.19) 

as A — >■ 00. It is convenient to study the asymptotic of this function using the representation (j3.14p . 
Notice that R{6 ± A) — > 1 as A — >■ 00. Besides, is straightforward to check that in this limit we have 

P'''''\Xe^,X'\Ye^,Y'\Ze^,Z') ~ P''{Xe^\Ye^\Ze^)P''' {X'\Y'\Z') (3.20) 
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Therefore, we get the cluster factorization property for the functions fj^^ , namely 

ff{9l, . . . , 6^+1 +A,...,9n+A)^ faiSn+l + A, . . . , f?Ar + A)/f (^i, ....On), (3.21) 

as A — > +00. The cluster factorization property for the form factors of the descendant operators 
immediately follows from (|3.13p . This result is in agreement with those of [TB]. The form factors of the 
level (n, n) descendant operators factorizes into the form factors of the chiral level (n, 0) and antichiral 
level (0,fi) descendant operators. 

Let us discuss certain consequences which follow from the cluster factorization property. We consider 
the chiral elements h € An and h' G An- The corresponding descendant operators will be denoted by 
and Vjj'' . The scaling dimensions of these operators in the ultraviolet limit are given by A^' = + n and 
Ajj = Aa + n, where Aq is the scaling dimension of the corresponding exponential operator. The spins 
of these descendant operators are the following, = n and = —n. From the cluster factorization 
property we immediately conclude that the element h corresponds to the chiral level n descendant 
operators while the element h' corresponds to antichiral level n descendant operators. 

Let us consider the generic descendant operator related to the element g = hh' . The scaling 
dimension of this operator is given by A^'' = Aa+n+n' while its spin is the following, s'*'' — n—n. From 
the cluster factorization property we conclude that this operators is given by certain linear combination 
of the level (I, I) descendants where I < n and I < n. Therefore, that the cluster factorization property 
do not impose rigorous conditions on the descendant operators which contribute to a particular function 
fg 



3.4- Descendants counting 

In this subsection we shall prove that the number of independent solutions to the form factor axioms 
in each level subspace coincide with the number of the descendant operators in the corresponding level 
subspace in the Lagrangian formalism. First, let us prove the following theorem 

Theorem 1. For generic values of the parameter a the map (I3.12p from the algebra A<E)A into the space 
of the functions fl is a bijection. 

Let us prove that the different elements g correspond to the different functions /^. First, we consider 
the chiral element g = h E A. Let us prove the linear independence of the set of polynomials (|3.17p . 
Indeed, for large enough N the functions 

z„ = S,n{X) - i-irS,n{Y) + (c^-™ - (-1)"l.")5™(Z) (3.22) 

are functionally independent. Therefore, the linear independence of the set of polynomials P'*(X|y|Z) 
reduces to the evident linear independence of the monomials z^^ . . . z^^ . 

Now, let us consider the asymptotic of the function ^ as the parameter a — > — ioo. It is easy to 
check that 

i7r(4\/2a-b+l/6) , iTr , , 

e 50 +-J^^^{xi,...,xn) ^iaix,)...aixN),h)^P''iX\0\0). (3.23) 

a— > — ioo 

This expression defines a map from the algebra A to the algebra of polynomials in the variables Zi — 
P''-{X\0\0). This map is evidently invertible. The linear independence of the polynomials P^{X\Y\Z) 
was already proved. Consequently, different elements from the algebra A correspond to different functions 
Now we apply the deformation argument. Since the map from the space of elements h E A to the 
space of form factors is a bijection at one point in the parameter a and the form factors are analytic 
functions in this parameter this map is a bijections for nearly all values of the parameter a. Hence, 
we prove that the form factors with different values of h differs. Besides, one can easily prove that 
dim An = dim Tn ■ 

Now let us consider the generic element g = hh' G A® A. Taking into account the cluster factorization 
property we conclude that the form factors of the different elements differs only if chiral and antichiral 
components of these form factors differs. Hence, we proved Theorem [T] 

As an immediate consequence of the theorem we have 

Proposition 1. For generic values of the parameter a the dimension of the space of the operators V£ 
with g e An'SiAn is equal to the dimension of the corresponding subspace of the Fock space dim(J>j(g) J^^). 
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3. 5. Integrals of motion 

The BuUough-Dodd model possesses a set of commuting integrals of motion Is of odd integer spin s 
except the multipliers of 3, i.e. 

s = 6n±l n = 0,1,2... (3.24) 
The local integrals of motion are diagonalized by the asymptotic states and the corresponding eigenvalues 
are given by 

N 

/,|0i,...,0A.) =^e^^'|0i,...,M, (3.25) 

where the proper normalization of the integrals of motion is assumed. On the other hand, let us consider 
the chiral element a-m- This element produces a common factor in all terms in (j3.15p . provided that 
the following relations are satisfied 

€""==-1, e-'3'"' + e^™ = 1. (3.26) 

The first relation shows that the value of m need to be odd, while the second one shows that the allowed 
values m are those of p.24p . For the corresponding elements a^s the common factor in all terms in p.lSp 
is of the form 

P"-=(X+|X_|Xo) = Ss{X+) + Ss{X^) + Ss{X^) = Ss{X), (3.27) 
where X = X- U X^ U Xq. Consequently, we get 

N 

/r'«(0i,...,0^)= ^e^^"/,s(0i,...,0jv), (3.28) 

■m— 1 

for any g ^ A® A. Comparing (|3.25p with (|3.28p we conclude that the elements a_s correspond to the 
appropriately normalized spin s integrals of motion Is while the antichiral element as correspond to the 
integral of motion i-e. 

V^-^^{x) = [VS{x),Isi V^-^<^{x) = [V^{x),I^s]. (3.29) 



4. Recurrent relations and reflection property for descendant operators 

In this section we shall prove that the form factors of descendant operators satisfy the refiection relations. 
These relations establish connection between the exponential operators with different values of the 
parameter a. Namely, up to some a-dependent factor the following operators are supposed to be coincide 

Va{x) = R{a) ^Q,-a(x), V-a{x) - i?'(a) V^Q,^+a{x), (4.1) 

where we introduce the notation 

The functions R{a) and R'{a) are the refiection amplitudes [Hi. These reflection relations originates from 
the similar relations for exponential operators in the Liouville theory [12J. Indeed, the Bullough-Dodd 
model can be treated as two different perturbed Liouville theories with charges or depending on 
which exponent in (|2.4I) is taken as the perturbing operator. By using recursion relations one can prove 
that the J functions of the exponential operators possess the following reflection properties [T7] 

JN,a{xi, . . . , Xn) = JN,QL-a{xi, . ■ ■ ,XAr), JN-aixi, ■ . . ,Xn) = J-Q'^^+aixi, ■ ■ ■ ,Xn)- (4.3) 

We argue that the correspondence between the exponential operators (|4.ip can be extended to the whole 
space of descendant operators. More speciflcally, for any descendant operator of the exponential field 
there exist a descendant operator of the exponential field V^^ such that the form factors of these 
operators coincide. Here we introduce the notation w for the element of the finite group W generated 
by the elements wi and W2 such that 

Wia = Ql — a, W2a = —Q'l^ — a. (4.4) 

To prove the refiection property for descendant operators we need to introduce auxiliary construction. 
This construction is an alternative free field representation for form factors or, to be more precise for the 
functions ^ defined in p.l4p . 



7 



4-1. The stripped bosonization 

From (|3.14p it follows that each form factor is proportional to the function ^ up to uniform factor 
which depends on the number of particles only. For the exponential operators, i.e. g = I, the free 
field representation for the functions JN,a was proposed in [17 . Let us briefly recall the construction. 
Thereafter we extend this free field representation to the case of descendant operators. 

Let us consider a Heisenberg algebra with the generators c?^, n € Z, n ^ 0. These generators satisfy 
the following commutation relations 

Mm> = 0, Kt, rf,T] = mA^S^+nfi, (4.5) 
where the coefficients are given by 



TT{b 

— "^UJ ^ cus -— 1 UUb — CUb 

D ^ o 

Note that 



=4a;^^cos — (^cos— -cos gg J- (4-6) 



A- ^ At„ - i-irA+. (4.7) 

Let a be a central element of the algebra. Besides, let |l)a be the vacuum state satisfying the annihilation 
conditions for the positive generators d^\l)a = (n > 0) such that a\l)a — a|l)o- The Fock space 
generated by negative elements {n > 0) from the vacuum will be denoted by D^. Similarly, we 
introduce the Fock space which is generated by positive elements {n > 0) from the vacuum a(l|- 
This vacuum satisfies annihilation conditions for negative modes and a(l|a = a(l|a- The Fock spaces we 
introduced possess decompositions into the level subspaces I?^ — ®^o^a'n ^^^"^ — ®'n=n'^an- 
Let us consider the exponential operators 



A±(z) = exp^ ^z-", A"(z) =: \+ {zuo-^)\- {zu) : . (4.8) 

Here we use the notation : . . . : for the normal ordering procedure which is defined by the relations 

X^{z')X^{z') =: A±(z')A^(z') :, 
A-(z')A+(z) = A+(z)A-(z') = /(±c.)/(Jc.2) : A+(z)A-(z') :, ^^'^^ 

Let us define the current 

t{z)^e"P\+{z) + e-'''P\-{z) + h\^{z), (4.10) 



which looks much the same as Lukyanov's one (|33j). The functions JN,a which determine the form 
factors of exponential operators are given by the matrix elements of the currents t{xi)^ namely 

Jjv,a(a;i, ■ • -.xn) = a(l|i(a;i) ■ .■t(xN)\l)a- (4-11) 

In order to obtain the free field representation for the functions ^ we need an additional construction. 
Let us consider two representations of the algebra A in the Heisenberg algebra, tt/j and tt^, defined as 
follows 

7rfl(a_„) = 7r^(a-.) ^ (n > 0). (4.12) 

It is easy to get the following commutation relations: 

[7rK(a_„),A±(z)] = (±)"+iz"A±(z), 

[^L(a_„), A±(z)] = ~(T)"+^z-"A±(z), (4.13) 
[7rfl(a_,„),7rL(a_„)] = -m(A+)-i(l + (-l)'")5,„_„,o 

We shall also use the notation 

a{h\^a{l\^R{h), \h)a^^L{h)\l)a. (4.14) 
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Consider the chiral element h E A. Taking into account the commutation relations (|4.13p we easily 
obtain that the J functions of chiral and antichiral descendant operators are given by the following 
matrix elements 

JnA^I, ■ ■ ■,Xn) = a{h\ti{xi) . . .tN{xN)\l)a, J^.^ai^l,- ■ ■,Xn) = {l\tl{xi) . . . (x N)\h) a- (4.15) 

Now, let us consider the generic element g = hh. We denote the corresponding matrix element as follows 

J'n'A^Ii- ■ ■,Xn) = a{h\ti{xi), . . . ,tNixN)\h')a- (4-16) 

Notice that this matrix element do not coincide with the function J^\{X). However, these functions 
can be related to each other. Indeed, by pushing the element nji^h) and TTL{h) from the definition of 
the vectors a{h\ and \h')a in the expression (I4.16P in the opposite directions we get that this function is 
given by certain linear combination of the functions J^\, i.e. 

...xn) = Y.'^Ta(^^' • • ■ '^^)- (4.17) 

i 

If in the l.h.s. of this expression we consider the elements from the level subspaces h G An and h' e An 
then in the r.h.s the sum is taken over the J functions corresponding to the elements hi G Am and 
h[ € An', such that Ui < n and fii < n. 

The vectors a{h\ and \h')a will be referred to as the 'physical' ones since their matrix elements can 
be expressed as a linear combination of J functions. These vectors form a subspaces in the spaces 
2?f and Vj^. The 'physical' subspaces can be decomposed into the direct sum over the level subspaces 
pfl.phys ^ (goo^^pfl,phys pL.phys ^ (^oo^^pL.phys^ Evidently that 

dim V^^'P^^'' = dim Pf'.vf = dim J-„ . (4.18) 

Further we shall use an alternative definition of the 'physical' subspaces. Namely, these subspaces can 
be defined as kernels of the operators Z3„ defined as follows 

Dn^uj--''d+ +uj--^''d-, (n^O). (4.19) 

It is easy to check the consistency of these two definitions. Indeed, the operators D„ satisfy the following 
commutation relations 

[Dn,7rRih)]=0, [Dn,nL{h)]^0, [i?„,i?„]=0. (4.20) 

Besides, one can show that for any state a{v\ E T)^ such that a(w|I?-„, > 0, there exist the element 
h E A such that a{v\ = a(l|7rfl(/i). Similarly, for any state \v)a E T>^ such that Dn\v)a — 0, n > 0, there 
exist the element h E A such that \v)a — TTL{h)\l)a- 

4-2. Analytic properties of the Jf^ ^ functions 

The form factors of exponential operators are known to satisfy recurrence relations |17| . One can obtain 
the similar relations for the form factors of descendant operators. Indeed, the function J%j^i ^{z,X) is 
an analytic function in the variable z depending on the parameters X — {xi, . . . , xat}. This function has 
simple kinematical and bound state poles at points prescribed by the form factor axioms. The residues 
at these poles can be evaluated explicitly as shown in [17] . One can separate the kinematical and the 
bound state poles contribution from the regular part and obtain 

N 

J%+U^,X) = f°°Xai^,X) + ^ ^l^Kn{X)J%_,^SXn) + 

71—1 

N 2 N _2 

+ E -^^B+{X)J^^JXnU^,Xn) - E /_"". B-iX)J^^JXnU.-\Xn), (4.21) 
n—1 n—1 
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Here we introduce the notation Xn = \ {x,i} and the functions Kn{X) and B^{X) are given by 

""■ ""• (4^22) 

The derivation of the expansion (|4.2ip repeats the main features of those one for the exponential operators 
and we omit it here. The regular part is given by the function j'j^^ ^{z, X). This function is regular 
everywhere except the points z — and z — oo. Since the sum over residues is of the order 0(z~^) as 
z ^ oo the asymptotic behavior of the function as a function of z is governed by the regular part 
j(oo)9. Further we shall also use another expansion of the function in the vicinity of z = which is 
given by 

N _i 

N _i _2 ^-12 

- E i"" \ , B+iX)Jl,Jxr.u^^X^) + J2 \ ^ B-{X)J%^S^^u^~\X^). (4.23) 

n=l " n=l " 

In this case the asymptotic behavior is governed by the function ^{z,X). From expressions (j4.2ip 

and (j4.23p we easily obtain that the regular parts of these expansions are related as follows 

- J^N^Aai^^X) = D%,aiX), (4.24) 

where the function Z?^ ai-^) given by 

N N 



D%JX) = ^nU;'B+iX)Jl,JxnU;; X)~Y. ^"^"'^n iX)Jl,JxnU;-';Xn) + 

N 

+ J2Kn{X)J^^_,jXn). (4.25) 



n— 1 n— 1 

AT 



The expansions (|4.2ip and (I4.23P allows us to calculate form factors recursively. The recursion relations 
for exponential operators, i.e. g = 1, immediately follows from these expansion since the regular parts 
and J*^""^ can be easily calculated [T7]- In Section [S] we use these expansions to obtain recurrence 
relations for the chiral level 2 descendant operators. 



4-3. Reflection relations for descendant operators 

In this subsection we shall prove that the form factors of descendant operators satisfy the reflection 
relations. More precisely, we prove this statement for the matrix elements (|4.16p . However, from (|4.17p 
it follows that if the reflection relations hold for J functions these relations hold for J functions as well. 

Theorem 2. For generic values of the parameter a there exist a representation ra{w) of the group W 
on the algebra A such that for any elements h,h' Cz A the following relation holds 

-JnA^I^-'-^Xn) - JN.wa [Xi,...,Xn)- (4.2b) 

The proofs of this theorem is very similar to those of [10 . The main idea is that all form factors 
can be obtained as a coefficients of the large rapidity expansion of the form factors of the exponential 
operators [TH], |19j . 

First, let us introduce auxiliary current 

s{z)=:X+{zio-'^^)\-{zio'^^):, 
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such that 

siz)t{x) = t{x)siz) = ffluj'^Afflco-'/^) : s{z)t{x) : . 

It is easy to check that the matrix element of some number of the currents t{xi) and s(yi) satisfy the 
reflections relations, namely 

a(l|t(Xl) . . .t{XK)s{y\) ■ . .S(2/L)|l)a = wa{AAx.{) . . .t{XK)s{V\) ■ • • | l)t«a- 

Indeed, the matrix elements of the currents t(xi) satisfy the reflection property while an insertion of the 
currents s(yi) in matrix elements results in the a-independent overall factor. 

Now let us prove the reflection property. First, we prove that the product of the currents t{xi) and 
s(j/i) acting on any highest weight vector spans the whole Fock module corresponding to this vector. 
Consider the expansion 

oo 

a(l|t(^f^2;) . . . t{i-^ z)s(r]{^ z) . . . s{rj-^z) = ^ 2;"" ^(n; ^1, . . . , ^p; 771, . . . , ?7,|. 

n=0 

Hereinafter we shall use the notation S = {£,1 , ■ ■ ■ , £,p) and H = (771, . . . , 77^). We propose that for generic 
values of the parameter a and large enough values of p and q it is possible to choose a set S^*^ and H^'^ 
with i = 1,. . .dim(J'®^)„ such that the vectors a(?i; S^*), H^'^j form a basis in the level subspace V^^. 
Let us prove this proposition in the limit a — >■ —ioo. In this limit we have 

a{im'^) ■ ..t{^;'z)s{T,^'z) . . . Sin^h) = 

= . . . • • -^VMM ■■ A+(Cr'-2) • • • X+{^;'z)s{v^'z) . . . s{ri-'z) : . 

Here -F(^i, . . . ,^p,rii, . . . ,r]q) is some irrelevant for our purpose function and the normal ordered product 
in the r.h.s. is of the form 



where 

1=1 j=i j=i 

Let us consider the expansion 

00 

: A+(Cr'^) • • • X+{^-'z)s{r,^'z) . . . s{rj;'z) := ^ " (-)(«;S;H|. 

n=0 

Then we get 

n s 

(_)(n;S;H| = a(l|E E '^m.-nsYliHuAm + '^n.dnj 

5=1 ni,...5ns>0 i=l 
ni + ...+ns=n 

with some nonzero constants C„j...„^. Notice that all possible products of enter the r.h.s. 

For large enough values of p and q the functions with 1 < i < n are functionally independent 
and can be considered as an independent variables. Besides, the monomials . . . are linearly 
independent. Consequently, for any set of the numbers A^^---'^"^ with s = 1, . . . ,n and ni, . . . ,ns > 0, 
m + . . . + rig = n, we have 

S (Ti,...,(Tti 

ni ,. . .n,5 

for some values of ^. Hence, the vector generated by the numbers A"^^-;-^'^ is not orthogonal to 

at least one vector generated by the product of k^. It means that there is no vector in the dim(J^®^)„- 
dimensional space orthogonal to all vectors generated by the product of for any values of 2 and H. 
Consequently, the vectors (_)(n;E!;H| for some values S^*) and H^'^ with i = 1, . . . dim(J^®^)„ form a 



11 



basis in the level n subspace of the Fock module. The deformation argument proves the same statement 
for the generic values of the parameter a. Hence we prove that the product of the currents t(xi) and 
s{yi) acting on any highest weight vector spans the whole Fock module corresponding to this vector. 

Now let us prove that the reflections act on the space of the matrix elements with respect to the whole 
Fock module. Let a{n; i\ — a{n; 5*^*-'; H'*-* |, i = 1, . . . dim(J^®^)„ be a basis vectors in the level n subspace 
of the right Fock module. Similarly, let the states \n;j)a be a basis vectors in the level n subspace of the 
left Fock module. As follows from the reflection properties for exponential operators, we have 

a{n;i\t{xi) . . .t{xN)\n;j)a = wa{n;i\t{xi) . . ■t{xN)\n; j)wa- 

Let us prove that the restriction of the whole Fock module to the 'physical' one preserves the reflection 
property. Let the vectors a(l|7r_R,(/ia.n,/^) — ainTi^l = '^i'^i{'^)a{n;i\ form a basis in the D^'^^^^. 
Similarly, let T^hih'g^ n v)\^)a — \n; v)a = X]j j)a be a basis in the T>^'^^^^^. We deflne the 'physical' 

states by means of the operator D,„ introduced in (|4.19l) . Notice that the auxiliary current s{z) is given 
by the following expression 

From this relation we easily get 

a{l\t{xi) . ..t{xM)s{yi) . . .s{yN)Djnt{x[) . . .t{x'M>)s{y[) . ■.s{yN')\l)a = 

wa{Mt{xi) . . . t{xM)s{yi) . . . s{yN)Dmt{x[) . . . t{x'M')s{y'i) ■ ■ . s{yN')\l)wa- 

We have 

= ainTlJ-lD-mln - m;j)a = ^^v^ {a)a{n\i\D_ra\n - m;j)a = ^v'j'{a)^a{n;i\D_m\n - m;j)nia- 

i i 

Therefore, 

(a)«)a(n-; j|£'-m =0. 

i 

From this equation we immediately conclude that there exist an element /i™^ „^ such that 

wa{l\T^R{hZa,n,,,) = ^V^ {a)wa{n;i\. 

i 

Similarly, there exist an element h'^^ „ ^ such that 

'^L{h'^a,n.,t,)\^)wa ^^v"^ {a)\n] j) ^a- 
3 

Finally, we get 

a{l\l^R{ha,n,tj)t{xi) . . .t{xN)TTL{h'a^fi^^)\l) a = a{n^\t{xi) . . . t{x N)\n, 1^) a = 

= ^'"i ("■)a{'>^'' '>'H^i) ■ ■ ■t{xM)\n\j)a = ^uf(a)uj(a)^a(n;i|t(a;i) . . .t{xM)\n\ j)wa = 

i,j hi 

«.a(lki?(/l™<,_„^^)i(xi) ...i(xjv)7rL(/l'™„ jj_^)|l)u,a- 

We obtain the map ra{w) on the subspaces p^^phys ^nd Vj^'^^y^, namely 

raiw)ia{l\TTRiha,n,t,)) = (1 1 TTfl (/l^„^„^^) , (w) (tTl (/i;,^^ 1 1) = L{h'ZMrn.u)\'^) wa- 

Comparing with the property we conclude that 
This proves Theorem [51 
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5. Explicit example: level (2, 0) descendants 

In this section we obtain recurrence relations for the form factors of descendant operators at the chiral 
level (2,0). Taking into account expansion (I4.2ip it is sufficient to find the regular part of the function 
Jf^ ^{z,X). This means that we need to consider the asymptotic of this function as z — oo. 

The level (2,0) subspace is spanned on the elements a_2 and a'^_^ acting on the highest weight 
vector a(l|- For these elements we get the following expansion 

J^+.^Jz, X) = &^1%{X) ■ + Citi(^) • ^ + ^So(^) + Oiz-'l (5.1) 
where the coefhcients are given by 

Q,a,i (^) = (2V3sm^p + /i)J^;; + ^— ^Ji,a • LnAX), 

C^a',2"^(^) = Jl^a ■ JnA^)^ 

^(oo)a_2/v^ iV^jh'^ - 1) J a_i n/3(/l^-l) . . ^^ 

^N,a,l \^)- 2 "^l.^'^JV.a "I 2 ' 

and 

L^Ax) = E h*^o,i#i.^#i-)^^r^Y: X. - E X, - ^ 5: X 

i++i-+io=i iei+ jei- ^ fcG/o 

X n /(f;")/(5-0/(S")/(S-') n f{^y 

iei+.jei-Meio J k k {p<q)eio * 

For the recurrence relations we also need the coefficients cj^°pJ^oiX) in the expansions (|5.ip . It is rather 
complicated to find them directly. However, it is not necessary since these coefficients are related with 
the leading coefficients in another limit z 0. The expansion of the ^ function in the vicinity of this 
point is given by 

Jl,^,Jz, X) = C'-^]1,{X) + d^AX) ■ z + d^]l,{X) ■ z' + 0{z'). 
By using factorization property (j3.2ip we immediately get the coefficients C^*']^q(X), namely 

As immediately follows from (|4.24p the required coefficients are given by 

Moo)g _ Mo)g n 

Consequently, expansion (|4.23p becomes a recurrence relation for the form factors of descendant operators 
at the level (2,0) where the regular part is of the form 

Jl,°il.(-, X) = &^ll{X) ■ z? + d^-^'AX) ■ z + dZ,{X). (5.2) 

Now let us find such combinations of the elements a_2 and which are invariant with respect to the 
reflections. It is easy to check that the required combinations are given by 



Q-2 
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(/i2-i)(y3 + 2/icos7r(p+l/6))V" ^ Ji,a ""^ 
Indeed, the element a-i corresponds to the integral of motion (I3.28p . consequently 



2 
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and the functions Jn.a are known to be invariant with respect to reflections. 

Now let us consider consider the element ha . We can easily obtain recurrence relation for this 
element. This relation is given by the expansion (j4.23l) where the regular part is given by the appropriate 
linear combination of the functions (j5.2l) , namely 

X) = -z^Jt: {X) + Jl,. . 42 (x). (5.3) 

This relation together with the initial conditions 

JqI = -^i.a (2;) = 0> Ji,a = 2cos7rp + 
defines the set of functions corresponding to the element ha uniquely. The invariance of this set of 



functions with respect to the reflections generated by the group W simply follows from (|5.3p . Indeed, 
one can easily check that the parameter a always enters the recurrence relations by means of the function 
cos np which is invariant with respect to reflections. 

The recurrence relations allows us to calculate the first few J functions at the level (2,0) easily, for 
example 

jhf^ , , r,-Ji,a-(rt + {Ji,a - h^ + h)ala2(^3 - Ji,a ' crj^l - /i(4 - 5h^ + h^)alal 

Y{i<j{^i+XiXj+xj) 

where we introduce the set of elementary symmetric polynomials cri(xi, . . . , xjq) defined by the generating 
function 

N N 

Y[{x + Xi) = ^a;^"Vi(xi, . . .,xn)- 

In principle, one can obtain recurrence relations for the form factors of descendant operators at each 
level. However, an explicit expression of these relations is rather cumbersome and not too useful. The 
recurrence relations provides a convenient framework to study analytical properties of form factor. 



6. Conclusion 

In this paper we use the free field approach to construct the space of solutions to the form factors axioms. 
We prove that there exist a bijection between this space and the Fock space of descendant operators 
for a generic values of the parameter a. The proposed free field representation provides a convenient 
framework to study analytical properties of the form factors. However, the problem of identification 
between the form factors and particular descendant operators is not solved. The cluster factorization 
property and the correct asymptotic behavior do not solve the identification problem except for the case 
of exponential operators. The reflection relations between the form factors might be a step forward to the 
solution of this problem. The further speculations about the identification problem involve the resonant 
identities between the operators PPI. 

The quantum group restrictions of the Bullough-Dodd model are known to describe certain classes of 
perturbed minimal models. It is natural to conjecture that the lightest breather form factors for any local 
operators in the perturbed minimal models can be analytically continued from the Bullough-Dodd case. 
We checked this conjecture for the form factors of primary operators in the Ising model in a magnetic 
field ^17j . Therefore, the proposed construction provides a convenient framework to study form factors 
of local operators in certain classes of the perturbed minimal models. 
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